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ABSTRACT 



We combine new Herschel/SPIRE sub-millimeter observations with existing multiwavelength data to investigate the dust scaling 
relations of the Herschel Reference Survey, a magnitude-, volume-limited sample of ~300 nearby galaxies in different environments. 
We show that the dust-to-stellar mass ratio anti-correlates with stellar mass, stellar mass surface density and NUV - r colour across 
the whole range of parameters covered by our sample. Moreover, the dust-to-stellar mass ratio decreases significantly when moving 
from late- to early-type galaxies. These scaling relations are similar to those observed for the Hi gas-fraction, supporting the idea 
that the cold dust is tightly coupled to the cold atomic gas component in the interstellar medium. We also find a weak increase of 
the dust-to-Hi mass ratio with stellar mass and colour but no trend is seen with stellar mass surface density. By comparing galaxies 
in different environments we show that, although these scaling relations are followed by both cluster and field galaxies, Hi-deficient 
systems have, at fixed stellar mass, stellar mass surface density and morphological type systematically lower dust-to-stellar mass and 
higher dust-to-Hi mass ratios than Hi-normal/field galaxies. This provides clear evidence that dust is removed from the star-forming 
disk of cluster galaxies but the effect of the environment is less strong than what is observed in the case of the Hi disk. Such effects 
naturally arise if the dust disk is less extended than the Hi and follows more closely the distribution of the molecular gas phase, i.e., 
if the dust-to-atomic gas ratio monotonically decreases with distance from the galactic center 

Key words. ISM; dust, extinction - Galaxies: evolution - Submillimeter; galaxies - Galaxies: clusters: individual: Virgo 



1. Introduction 

One of the main challenges for extragalactic astronomy is to 
understand how galaxies evolved from simple clouds of un- 
enriched gas into the complex eco-systems harboring stars, plan- 
ets, dust, metals and different gas phases we observe today. Of 
course, this transformation must be related to the star forma- 
tion cycle in galaxies. The gas cools, condensing into molecular 
clouds that collapse forming stars. Heavy elements are formed 
during nuclear reactions in the stellar cores and then, when the 
stars die, they are expelled into the inter-stellar medium (ISM) 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



and either mix with the gas phase or condense forming dust 
grains. However, we still know very little about the detailed 
physical processes regulating the star formation cycle in galaxies 
and how much it is affected by the internal properties of galax- 
ies (e.g., mass, dynamics and morphology) and the environment 
they inhabit. 

In the last decades, wide-area optical surveys of the local 
Universe and deep investigations of high-redshift systems have 
allowed an accurate reconstruction of the average star formation 
history (SFH) of g alaxies in different environments across th e 
Hubble time (e.g.. lAbazaiian et~al ] 120091: IScoviUe et al ] l2007h . 
They have not only revealed that the star formation rate (SFR) 
density of the Univer se had a peak around z ~1 and is now 
rapidly declining (e.g.. lLillv et"al]ll996l:lMadau et alJll998h . but 
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they have also shown that the SFH of galaxies is tightly linked 
to their total mass (i.e., massive systems formed the bulk of 
their s tellar populations ea rlier than dwarf galaxies; Co wie et al.l 
Il996t lOavazzi et anil996l) . and to the environment (i.e., at the 
current epoch, galaxies i n high density regions are less ac- 
tive than isolated system s; [Lewis et al. l l2002l:lG6mez et al l2003t 
iBoselli & Gavazzi l^2006'). The next step is to investigate why the 
conversion of gas into stars is mass and environment dependent 
and how this impacts on the metal enrichment of the ISM. 

A unique role in the star formation cycle of galaxies is played 
by dust grains. Since dust is formed out of the metals produced 
during stellar nucleosynthesis, its properties provide us with im- 
portant clues about the recent star formation activity of galaxies. 
Moreover, dust is also believed to act as catalyzer for the for- 
mation of molecular hydrogen (H2) and prevents its dissociation 
by the interstellar radiation field, thus helping to regulate the fu - 
ture star formation activity (e.g.. lHollenbach & Salpete3ll97lh . 
Finally, since dust can also be de stroyed, in particular if not em- 
bedded into the cold ISM (e.g.. IClemens et alll201(]|) . its pres- 
ence and abundance can yield important information about the 
properties of the other components of the ISM in galaxies. 

Alth ough dust properties h ave been studied for several 
decades dSavage & Mathisll 19791) . it has become gradually clear 
that, in order to obtain a proper census of dust in the Universe, 
it is necessary to incorporate extragalactic investigations of the 
far-infrared/sub-millimeter (submm) regime. It is in this part 
of the electromagnetic spectrum that the cold dust component 
dominating the dust mass budget in galaxies hides. Pio neers 
in this area have be en the IRAS ( Neugebauer et al J 19841) . ISO 
dKessler et al.ll 19961) and Spitzer (Werner et al. 2004) telescopes 
and the SCUBA can iera (iHoUand et al. 199 9). but now with the 
launch of Herschel (iPilbratt et alJ 20101) are we entering a new 
era for interstellar dust studies. Thanks to its unprecedented sen- 
sitivity in the wavelength range 200< A <600 fim and a much 
improved spatial resolution at 70< A <160 fim (with respect 
to Spitzer), Herschel is a unique instrument to unveil the role 
played by dust in the evolutionary history of galaxies. 

The first, natural step in this direction is to quantify how 
the dust content of galaxies varies with internal galaxy proper- 
ties such as, stellar mass, colour, surface density, gas content, 
etc. These scaling relations will provide initial clues on the role 
played by dust on the star formation cycle and strong constraints 
for chemical evolution models. For example, it is important to 
quantify the link between the dust and the cold gas component 
of the ISM and investigate whether or not they follow similar 
scaling relations. 

The main scaling relation investigated during the pre- 
Herschel era has been the evolution of the dust- to-g as ra- 
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with stellar mass and metallicity (e , -- 
iLisenfeld & Ferrara 1998; iPopescu et al.l I2OO2F 
20071; iGalametz et al.ii201 li) . Several studies have shown an in- 
crease of the dust-to-gas ratio as a function of stellar mass, 
mimicking what has already been observed in the case of the 
stella r mass vs. gas-phase metallicity relation jTremonti et al.l 
I2OO4I) . Although several investigations have quantified the vari- 
ation of the total far-infrared-to-optical/near-infrared luminosity 
(usually interpreted as a proxy for the ratio between obscured 
star formation rate and stellar mass) as a function of morpholog- 
ical type (e. g., 'de Jong et aL[ Il984t iDevereux & Hameedlll997t 
iBendo et al.l l2002; Popescu et al. 20021) . surprisingly verv little 



is known about the relation between the dust-to-stellar mass ratio 
and galaxy properties. Only very recently, da Cunha et al. (2010) 
have shown that the dust-to-stellar mass strongly correlates with 
specific star formation rate (i.e., the ratio between the current 



star formation rate and the stellar mass, SSFR), as predicted by 
chemical evolution models. However, it is still unknown whether 
or not morphology and stellar mass also regulate the amount of 
dust present in nearby galaxies. 

Even more limited is our knowledge of the role (if any) 
played by the environment on the dust content of galaxies. 
Whereas it is now well established that the atomic hydrogen 
(Hi) content of galaxies depends on the environment they in- 
habit, it is still debated whether or not the harsh cluster envi- 
ronment can also affbct the molecular hydrogen c omponent in 
the ISM (iBoselh et al.ll2002l; lFumagal li et al.l2009h . Thus, since 
the dust is supposed to be more directly linked to the molecu- 
lar than to the atomic gas phase, it is not clear if dust is removed 
from the star-forming disks of cluster galaxies. Before the launch 
of Herschel, no definitive evid ence of dust stripping in cluster 
galaxies had been found (e.g.. Tuffs et al. 2002; Popesc u et al.l 
i2002; Contursi et al. 2001; BoselU & Gavazzi 2006).' 

The situation has considerably improved in the last year, 
thanks to the Herschel Science Demonstration Phase. Herschel 
observations have not only allowed u s to directly witness 
dust^ strip ping in clusters of galaxies dCortese et al.l 1201 Oat 
IGomez et aLlbOlOt) . but they have also clearly demonstrated that 
the dust disk is significantly reduced in very Hi-deficient cluster 
galaxies, fo Uowiiig remarkably w ell the obs erved 'truncation' of 
the Hi disk dCortese et alj201 0b'). Moreover. IWalter et alj d2011h 
have found evidence of dust stripping also in the M81 triplet. 
Thus, while it is becoming clear that dust can really be perturbed 
by the environment, it is still uncertain whether these are just ex- 
treme cases or dust stripping is a common evolutionary phase for 
cluster galaxies. 

Dust scaling relations are the necessary step to solve this is- 
sue: only after the relation between internal galaxy properties 
and dust content has been properly quantified, will it be possi- 
ble to determine the role played by the environment. It is in fact 
mandatory to compare the properties of galaxies in different en- 
vironments at fixed stellar mass, morphological type and colour 
to exclude that our findings are just a result of secondary trends 
between galaxy structure, star formation activity and environ- 
ment. 

For all these reasons, in this paper we take advantage of 
Herschel observatio ns of the Herschel Reference Survey (HRS, 
iBoselli et al.l2010bl) to quantify how the dust-to-stellar mass and 
dust-to-Hi mass vary with internal galaxy properties and envi- 
ronment. The HRS is the ideal sample to carry out this analy- 
sis in the local Universe. Its selection criteria (magnitude- and 
volume-limited), size (~300 galaxies), sensitivity to the cold 
dust component (down to ~10"' M© at the distance of the Virgo 
cluster) and multiwavelength coverage, make this sample ideal 
to investigate dust properties in the local Universe, thus provid- 
ing strong constraints for theoretical models and a reference for 
high-redshift studies. 

This paper is organized as follows. In § 2 we describe the 
Herschel observations and data reduction, the technique used to 
estimate the dust mass and the ancillary data used to trace the 
other components of the ISM. In § 3 and § 4 we discuss the rela- 
tion of the dust-to-stellar mass and dust-to-gas mass with galaxy 
properties and environment while in § 5 we discuss the impli- 
cations of our results for chemical evolution models and for en- 
vironmental studies. In Appendix we investigate possible biases 
affecting our dust mass estimates. 
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2. The data 

The HRS consists of a volume-hmited sample (i.e., 15< D <25 
Mpc) including late-type galaxies (Sa and later) with 2MASS 
(Skrutskie et al. 2006) K-band magnitude K^^, < 12 mag and 
early-type galaxies (SOa and earlier) with Y^stot ^ 8.7 mag. 
Additional selection criteria are high galactic l atitude (b > +55° ) 
and low Galactic extinction (A^ < 0.2 mag, Sch legel et al.l 19981) . 
to minimize Galactic cirrus contamination. The total sample 
consists of 322 galaxies (260 late- and 62 early-type g alax- 
ie£l). As extensively discussed in Bose lli et al.1 (1201 Obi) , this 
sample is not only representative of the local Universe but also 
spans different density regimes (i.e., from isolated systems to 
the center of the Virgo cluster) and so it is ideal for environmen- 
tal studies (see also Hughes & Cortese.2009; Cortese & Hu^hesl 
2009; Boselli et al jj2009HCortese etal.ll201lh . As discussed in 
Boselli et al.i (i201()b l). we fixed the distances for galaxies belong- 



ing to the Virgo cluster (i.e., 23 Mpc for the Virgo B cloud and 
17 Mpc for all the other clouds; iGavazzi et al.1 ^1999). while for 
the rest of the sample distances have been estimated from their 
recessional velocities assuming a Hubble constant Hq =70 km 
s"' Mpc"'. 

2.1. Herschel/SPIRE observations, data reduction and flux 
extraction 



Herschel/SPIRE (Grifl in et al.ll2010h observations of galaxies in 
the HRS were originally planned as part of the Guaranteed Time 
of the SPIRE and PACS consortia. However, after the selection 
of the Herschel Open Time Key Projects, it emerged that 83 HRS 
galaxies were also targeted, to a si milar de pt h, by the Herschel 
Virgo Cluster Survey (HeViCS, iDavies et am2010l) . a key pro- 
gram focused on the study o f ~60 deg^ in the Virg o cluster using 
PACS/SPIRE dGriffin et al.ll2010t iPoghtsch et al.ll20Toh parallel 
mode. In order to reduce as much as possible the duplication 
of Herschel observations, the local galaxies group of the SPIRE 
science team (SAG2) and the HeViCS team have agreed to share 
the data obtained for the galaxies in common between the two 
programs. Thus, the Herschel/SPIRE data presented in this work 
have been obtained as part of these two Herschel Key Projects. 

The 239 galaxies outside the HeViCS footprint, plus 4 galax- 
ies in HeViCS targeted by the HRS during the Herschel Science 
Demonstration Phase, have been observed using the SPIRE 
scan-map mode with a nominal scan speed of 30 "/sec. The size 
of each map was chosen on the basis of the optical extent of the 
target. Galaxies with optical diameters smaller than ~3' were ob- 
served using the small scan-map mode providing homogeneous 
coverage over a circular area of ~5' diameter For larger galax- 
ies, the large scan-map mode was used, with each map covering 
at least an area with diameter 1.5 times the optical diameter of 
the target. The typical sizes of these maps are 8' x 8', 12' x 12' 
and 16' x 16'. The depth of each map was decided according 
to the morphological type of the target: 3 and 8 pairs of cross- 
linked scan-map observations for late- and early-type galaxies. 



' Wi th respect to the original selection described in iBoselli et alj 
( 1201 Obh . we removed HRS228 because the redshift reported in NED 
turned out to be incorrect. In addition, after visual inspection and com- 
parison with the literature, we revised the NED-b ased morphologi- 
cal type classification given in IBoselli et alj (20 10b) for six galaxies 
in our sample. In detail, NGC4179 moved from Sb to SO, NGC4438 
from so/a to Sb, NGC4457 from SO/a to Sb, VCC1549 from Sb 
to dE, NGC4691 from SO/ a to Sa and NGC5701 from SO/a to Sa 
( Sandage & Tammann 1981; iBinggeli et allll985l ; IGavazzi et aLll2003l ; 
PCortese & Hughes,.2009i) . 



respectively. This resulted into a pixel-by-pixel rms at 250, 350 
and 500 jum of ~7, 8, 9 mJy/beam for late-type galaxies and ~6, 

6 and 7 mJy/beam for early-type systems. For additional details 
see Ciesla et al. (in prep.). 

The remaining 79 HRS galaxies targeted by the HeViCS sur- 
vey have been observed by Herschel using the SPIRE/PACS 
parallel scan-map mode with a fast scan speed of 60 "/sec. 
Eight cross-scans fields of ~ 4 x 4 deg^ each across the Virgo 
cluster were obser v ed co vering a total area of ~60 deg^ (see 
iDavies eta D 120101 1201 2' for additional details). The typical 
pixel-by-pixel rms measured around the HRS galaxies are ~7, 

7 and 8 mJy/beam at 250, 350 and 500 fim. 

Both datasets were reduced using the same data reduction 
technique. Extensive details of our data reduction procedure can 
be found in Smith et al. (in prep.). In the following we briefly 
summarize the main steps. The SPIRE photometer data were 
processed up to Level- 1 with a script adapted from the official 
SPIRE pipeline. The main differences with respect to the stan- 
dard pipeline are that we use a different deglitching method and 
did not run the default temperature drift correction and the resid- 
ual baseline subtraction. Instead, we used the BriGAdE method 
(Smith et al. in prep.) to remove the temperature drift and bring 
all bolometers to the same level (equivalent to baseline removal). 
This method provides a better baseline subtraction, in particular 
when strong temperature variations are present. All scans were 
then combined into a map using the naive map maker included 
in the standard pipeline. The pixel sizes of each map are 6, 8 
and 12 arcsec and the corresponding full widths half maximum 
(FWHM) 18.2", 24.5", 36.0" at 250, 350, 500 jum, respectively. 

Total integrated flux densities were obtained as extensively 
described in Ciesla et al. (in prep.). Briefly, for extended source^ 
elliptical apertures with a typical major axis 1.4 times the optical 
diameter of the galaxies were used. In some cases (e.g., early- 
type galaxies) the size of the aperture was significantly reduced 
in order to include all the emission from the galaxy and/or avoid 
as much as possible contamination from other sources or fore- 
ground cirrus emission. Background was estimated using circu- 
lar annuli around the target (Ciesl a et al ., in prep.). Errors were 
estimated following IBoselli et alj (l2003l) . For point-like sources 
(8 out of 322 galaxies at 500//m), flux densities were extracted 
directly from the time-line data using the method developed by 
Bendo et al. (in prep.). This turns out to provide a more reliable 
flux density estimate than the aperture photometry technique. 
Errors for point sources come directly from the time-line fitting 
technique. The typical errors on the flux estimates for all HRS 
galaxies are -6%, 8% and 11% at 250, 350, 500 yum, respec- 
tively. These do not include the SPIRE calibration error of ~7%. 
In case of non-detections, upper-limits have been determined as- 
suming a 3cr signal over a circular aperture of radius 0.3, 0.8 and 
1 .4 the optical radius for E, SO and spirals, respectively (Ciesla 
et al., in prep.). In the following, we will work with Relative 
Spectral Responsivity Function (RSRF)-weighted flux density 
measurements. Since the SPIRE pipeline produces monochro- 
matic flux densities assuming that the source is point-Uke, we 
had to apply correction factors to convert monochromatic flux 
densities into RSRF-weighted flux densities. Thus, we divided 
monochromatic flux densities obtained from the maps by 1 .0 1 1 3, 
1.0087 and 1.0065, at 250, 350 and 500yum, respectively (see 



' We define a source as 'extended' if a Gaussian fit to the time-line 
data provides a FWHM larger than 20", 29", 37" at 250, 350, 500 yum, 
respectively. 
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iHerschel Space Observatorvll201 f L and the SPIRE Photometer 
Cookboolfl). 

In total, 276 out of the 322 HRS galaxies (-86%) have been 
detected in all the three SPIRE bands, 38 galaxies are non- 
detections in all bands, 3 have been detected at 250yum only and 
5 galaxies at 250 and 350jum only. 

2.2. Dust mass estimate 

The main aim of this paper is to investigate the relation between 
dust content and internal and environmental properties of galax- 
ies. Thus, an accurate estimate of dust masses is mandatory. 
Ideally, the best method to quantify the amount of dust in galax- 
ies would be via spectral energy distribution (SED) fitting of data 
across the whole ultraviolet to far-infrared/submm regime us- 
ing up-to-date dust models and energy balance/radiative transfer 
techri i ques (e.g.,fe ubko et al. 2004; Draine & Li 2007; Bianch 



2008 



20IC; 



Galliano et al. 2008; jda Cunha et alj i2008i; iBaes et al 
Compiegne et al. assuming simple sin- 

gle/double modified black-body emission (e.g., iDunne et al.l 
I2OOOI; IVlahakis et aP '2005). Unfortunately, no homogeneous 
dataset for all the H RS galaxies is currently available for A <200 
fim (see'Boselli et al. 2010b), making an accurate SED fitting of 
only the three SPIRE flux densities unfeasible in the case of de- 
tailed dust models and perhaps not ideal even in the cas e of a sin- 
gle modified black-body emission ("S hettv et alJ2009allbh . This is 
mainly because, with only the SPIRE points, no constraints on 
the shape of the SED below 250 yum are available and no de- 
tailed error analysis can be performed. Thus, in order to avoid 
introducing systematic biases due to the incomplete coverage 
of the HRS in the mid- and far-infrared wavelength regime, we 
derive empirical recipes to estimate dust masses from the three 
SPIRE flux densities only. This approach is similar to what has 
been done in the past to determine stellar mas ses from optical 
and near-infrared l uminosities and colours ("e.g. jBell et al.l2003t 
IZibetti et alj|2009l) . In recent years, these empirical calibrations 
have become the ideal tool to roughly quantify the stellar con- 
tent of galaxies when accurate SED fitting is not possible. Thus, 
the technique here adopted must not be considered equivalent to 
a proper SED fitting, but it just provides us with a way to have a 
first view on the dust properties of the HRS sample. 

As a first step, we test the feasibility of this technique by 
using two different approaches. Firstly, we assume that the dust 
SED is well approximated by a simple single modified black- 
body radiation, then we will compar e our results with the pre- 
dictions of the models developed by Draine & Lil (|2007[) . A dis- 
cussion on the possible biases affecting our dust mass estimates 
is also presented in Appendix A. 

2.2.1. Modified black-body 

In this case, the dust mass is given bv (lHildebrandl[T983l) 



Mdust = 



KyBy(T) 



(1) 



10 



"T 1 1 1 1 — I I I 




• Black-Body 

• Draine & Li 



J I I I I I 



1015 



10'6 



Mdust/(f35oD') [kg W-1 Hz] 



Fig. 1. The relation between the ratio Mdust/iFisoD ) and the 
250yL(m-to-500yL(m flux density ratio for a modified black-body 
SED with /3=2 (black, see Table Bl) and for the Draine & Li 
models (red). 



Planck function 



By(T) = 



2hv^ 



1 



„hv/KT _ I 



(2) 



where h is the Planck constant, c the speed of light and T is the 
dust temperature. We can use EqUto determine how the 250-to- 
500//m flux density ratio varies with temperature: 



/2 



./250 _ / > 
/5OO 



' V5OO 



ghv25o/KT _ I 



(3) 



and thus, by combining Eq. [T] and [3] we can find a functional 
form that relates, for example, the ratio MdustlihsoD^) to the 
250-to-500 yum ratio. We tested this by varying the dust temper- 
ature in the range 5-55 K, a ssuming a dust mass absorption co- 
efficient of 0.192 at 350//m ( lDraindl2003l) and yS=2. We choose 
P-l for two reasons. Firstly, a r ecent analysis of V irgo galax- 
ies (some of which in our sample: Ibavies et ani2012h has shown 
that their dust SEDs are well fitted by a single modified black- 
body with y6=2. Secondly, this value is very c lose to the one 
assumed in the models by iDraine im (l2007b . thus simpUfy- 
ing the c omparison of our res ults with more refined dust models 
(see also lMagrini et al.ll201 ih . The result of this test is shown in 
Fig.[Tl where the black line has been obtained by convolving the 
black-body SED with the SPIRE Relative Spectral Responsivity 
Function (RSRF) for extended sources. Unsurprisingly, the two 
quantities are strongly correlated and it is easy to fit a polynomial 
function to the data, suggesting that this method could be easily 
used to estimate dust masses when SED fitting is not possible. 



where fy is the flux density emitted at the frequency v, Ky is 
the dust mass absorption coefficient [in the following we will 
assume Ky - /C35o(v/v35o)''], D is the distance and By(T) is the 



^ http://herschel.esac.esa.int/Docs/SPIRE/html/spire_om.html 
* http://herschel.esac.esa.int/twiki/pub/Public/SpireCalibrationWeb/ 
SPIREPhotometryCookbook.may 20 1 1 . 1 .pdf 



2.2.2. Tine Draine & Li (2007) model 

In order to check the validity of our method, we tested it by 
using a more detailed descriptio n for the du st SED as given by 
the dust models developed bv .Draine & Lil (1200 7 iFl According 

' http://www.astro.princeton.edu/~draine/dust/irem.html 



4 



Cortese et al.: The dust scaling relations of the Herschel Reference Survey 



to these models the dust mass in a galaxy can be related to the 
flux density emitted in the far-infrared/sub mm by the following 
equation 

^ Mh ' Jy 

where MdustlMn is the mass of dust per nucleon (given by the 
dust model assumed), nin is the mass of the hydrogen nucleon 
and /'v is the dust emissivity per H nucleon. As discussed by 
iDraine & Lil (I2007) . the dust emissivity strongly depends on the 
properties of the starlight intensities responsible for the dust 
heating. In particular, they show that the total dust emissivity 
can be expressed as 

jv = (1 - y)jv(Umm) + yjviUmin, U,nax, «) (5) 

where t/„„„ and Umax are the lowest and highest possible in- 
tensities for the interstellar radiation field, (l-y) is the frac- 
tion of dust mass exposed to starlight intensity f/„„„ and a 
is the index characterizing the distribution of starlight intensi- 
ties. They find that the SEDs of galaxies in the SINGS sam- 
ple appear to be satisfactorily reproduced with a fixed a -2 
and Umax =10^, so we use these parameters in the rest of 
our analysifl By using the models for a Milky-Way dust (i.e., 
from MW3. 1_00 to MW3. 1_60, seven models in total), we com- 
puted jy from Eq|5]by varying y between and 0.3 (step 0.01) 
and assuming all the possible values f or t/,„,„ (from 0.1 t o 
25, twenty-three in t otal), see als o Munoz -Mateos et a n (I2009h : 
IWiebe et alj (|2009|) : IWilkner et ah (i2009i) . Then, we estimated 
the 250-to-500 //m flux density ratio and, using Eq.[4l the ratio 
MdustKf35oD^) assuming MdustlMn =0.01 dPraine & Lill2007h 
and rriH =1.67x10"^^ kg. The red points in Fig. [T] show the 
relation between M^ust / (hsoD^) and the 250yum-to-500yum flux 
density ratio so obtained. We recover a relation extremely simi- 
lar to the one found for the modified black-body case (see also 
iMagrini et alJ201 H) . There is almost no systematic offset (~0.01 
dex) between the two methods and the typical scatter between 
dust mass estimates using the two approaches is ~0.02 dex (see 
also Appendix A). We have also compared the results obtained 
from this method with the dust masses estimated via SED fitting 
for those galaxies in the HeViCS footprint for which PACS 100 
and 160//m are available. We find that the two techniques pro- 
vide consistent results with no systematic effects and a typical 
scatter of ~0. 10 dex (see also Appendix A). 

In summary, we have developed a simple method that allows 
us to estimate dust masses from just the three SPIRE flux densi- 
ties. In the rest of the paper we will use the dust masses obtained 
from the 350/vm flux density and the 250-to-500 fim ratio as- 
suming a modified black-body with /3-2. For the non-detections, 
upper limits to the dust masses are computed by using the av- 
erage /og(/25o//50()) colour observed for the detections (0.75). 
Given the uncertainty in each flux density measurement, the typ- 
ical statistical error in the dust mass is ~0.2 dex (not including 
systematics due to assumptions on the dust properties). For the 
estimate of the dust mass, we also treated as non-detections those 
objects for which the submm emission is clearly aff ected by syn- 
chrotron emission: i.e., M84 and M87 (see also iBoselli et alj 
I2010al: I Baes et al] 12010) . The exact functional forms used for 
extended and point sources, as well as the one obtained for dif- 
ferent values of jS are presented in Appendix B. 



^ We also tested the effects of varying (/„„„ and we do not find sig- 
nificant differences. 



2.2.3. The effect of the dust emissivity on the dust mass 
estimate 

We have shown that our empirical method is consistent with a 
typical SED fitting technique, once the properties of the under- 
lying dust population are fixed. The iDraine et alj ( l2007l) model 
is consistent with j6=2 at these wavelengths, so it is not com- 
pletely surprising that we obtain similar results when com- 
pared with a modified black-body with /3-2. Unfortunately, 
the absolute value of the dust opacity and its depe ndence 
on fr equency are still quite uncertain (e.g., [ Finkbein er et al.l 
ITggg'-.'Dupac et al."2003':'Gordon et alj2010l:lParadis et alj2010l: 
Planck Collaboration et al. 201 1 ) and it is therefore crucial to in- 
vestigate how this assumption can affect our analysis. It is easier 
to address this issue by discussing separately a change in the 
absolute value of dust opacity and in the /3 parameter. A vari- 
ation in the absolute value of /C350 would just shift all our dust 
estimates up and down without affecting any correlation (and 
its significance) that we may find. In other words, this will just 
move systematically the y-axis in all our plots. Much more crit- 
ical would be a change in the value of /3, since this could not 
only create systematic offsets in the dust mass estimates but also 
alter the shape of the scaling relations by introducing secondary 
effects with the quantities we are interested in (e.g., stellar mass, 
colour and stellar mass surface density), thus undermining the 
reliability of our analysis. Therefore, in Appendix A, we have 
investigated how our analysis would be affected by using other 
values of namely, 1, 1.5 and 2.5. We find that, in addition to 
a systematic shift in the dust mass values (up to ~1 dex when 
moving from (3-2.5 to 1), a change of /? would also introduce 
some artificial trends, in particular with stellar mass and stellar 
mass surface density. In particular, moving from [3-2.5 to [3-1 
dust masses would gradually decrease, but galaxies with higher 
stellar masses and stellar mass surface densities would be more 
(~0.2 dex) affected than lower mass systems. Luckily, such dif- 
ferences between high and low stellar mass systems starts to be 
important only for [3 < 1 .5. For y6 > 1 .5 the scatter in the relations 
is <0.1 dex, i.e. comparable to, or lower than, the scatter in our 
method and, as we will see, significantly lower than the typical 
dynamical range covered by the scaling relations we investigate 
here. Moreover, it is important to note that the bluest submm 
colours in our sample (i.e., /250//500 >7.5) cannot be reproduced 
by a s ingle modified black-body with [3-\ (see also iBoselli et al.l 
l20T2l) . 

Thus, we can confidently proceed by estimating the dust 
masses using the technique developed above for a modified 
black-body with [3-2, and we will discuss throughout the text 
if and how much different values of [3 could influence our con- 
clusions. Given all the caveats discussed above, we remind the 
reader to be cautious of the absolute values of dust masses pre- 
sented here, and to mainly focus on the trends and diflFerences 
between the various samples discussed in the rest of the paper. 

2.3. Ultraviolet, optical and HI data 

The //er.sc/!e//SPIRE data have been combined with multi- 
wavelength observations covering the ultraviolet, optical and Hi 
regime. A detailed descript ion of these d a tasets and of the de- 
rived quantities is given in ICortese et al.l (1201 li) . Briefly, opti- 
cal broad-band photometry has been o btained from the Sloan 
Digital Sky Survey DR7 (SDSS-DR7, 'Abazajian et al.' '2009|) 
database. Galaxy Evolution Explorer (GALEX, Martin et al] 
■2005) near-ultraviolet (NUV; ^=2316 A: Ai=1069 A) images 
have been mainly obtained as part of two on-going GALEX 
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Fig. 2. The K-band luminosity distribution for the whole HRS 
(black) and for the sub-sample analyzed in this paper (blue). 



Guest Investigator programs (GI06-12, P.I. L. Cort ese and GI06- 
01, G ALEX Ulti-aviolet Virgo Cluster Survey, i Boselli et all 
1201 Ih . Additional frames have been obtained from the GALEX 
GR6 public release. The SDSS images were registered to the 
GALEX frames and convolved to the NUV resolution. Isophotal 
ellipses were then fit to each image, keeping fixed the center, el- 
lipticity and position angle (generally determined in the /-band). 
The sky background was determined within rectangular regions 
around the galaxy and subtracted from the images before per- 
forming the ellipse fitting. Asymptotic magnitudes have been 
determined from the growth curves obtain ed following the tech- 
nique described by iGil de P az et al.' (2007) and corrected for 
Galactic extinction assuming a Cardelli et al. (1989) extinction 
law. An extensive description of GALEX and SDSS measure- 
ments as well as all the fluxes will be presented in a future work 
(Cortese et al. in prep.). Stellar masses M» are determined from 
/-band luminosities L, using the ^ — i colour -dependent stellar 
mass-to-light rati o relation from 'Zib etti et a l.| 1*2009), assuming 
a lChabrieTl(f2003h initial mass function (IMF). We assume a typ- 
ical uncertainty of 0.15 dex in the stellar mass estimate. 

Atomic hydrogen masses have been estimated from Hi 21 
cm line emission data (mainly single-dish), available from 
the literatu r e (e.^ 
iKent et alJ 1200? 



Sorinsob et a 


J 


2005; 


Gavazzi et al. 




2003 



and die NASA/IPAC 

Extraealactic Database, NED). As described in ' Cortese et alj 
(1201 ih . we estimate the Hi deficiency parameter (Defni, defined 
as the difference, in logarithmic units, between the expected Hi 
mass for an isolated galaxy with the same morphological type 
an d optical diameter of th e target and the observed value, follow- 
ing 'Ha\Ties'&Giovaneni[i984) using coefficients that vary with 
morphological type. These coefficients have been calculated for 
the foUowing ty pes (see also Table 3 inlBoselli & Gavazzil2009l) : 
SOa and earUe r (IHavnes & Giovanellill984 . Sa-Sab, Sb, Sbc, Sc 
(Solanes etalJll996l) and Scd to later types dBoselU & Gavazzil 
[2009 ). It is important to remember that, for morphological types 
earlier than Sa, t he estimate of the exp ected Hi mass is highly un- 
certain (see also lCortese et alj|201 Ih . In the following, we will 
consider as 'Hi-deficient' galaxies those objects with Defni > 



0.5 (i.e., galaxies with 70% less hydrogen than isolated systems 
with the same diameter and morphological type). 

In order to avoid any possible bias in the comparison among 
the various scaling relation, in this paper we focus our atten- 
tion only on those galaxies for which Herschel as well as Hi, 
NUV and SDSS observations are currently available: 282 galax- 
ies (~87% of the whole HRS, namely 234 late- and 48 early-type 
galaxies). In Fig.|2l we compare the K-band luminosity distribu- 
tion of the whole HRS (black) and the one of the sub-sample here 
investigated (blue). It is clear that our sub-sample provides a fair 
representation of the HRS and therefore the scaling relations in- 
vestigated in the rest of the paper should be representative of the 
local galaxy population as a whole. 

3. The dust-to-stellar mass ratio 

As a first step in order to quantify how the dust content varies 
as a function of integrated galaxy properties, we plot in Figl3] 
the dust-to-stellar mass ratio as a function of stellar mass M* 
(upper-left panel), observed NUV - r colour (upper-right), con- 
centration index in r-band [C3i(r), defined as the ratio between 
the radii containing 75% and 25% of the total r-band lighQ] and 
stellar mass surface density /v, [i.e., Mt/(2nR^Q .), where Rsoj 
is the radius containing 50% of the total /-band light]. Circles 
and triangles indicate detections and upper-limits, respectively. 
The dust-to-stellar mass ratio strongly anti-correlates with M,, 
NUV - r (a proxy for sp ecific star formation rate, SSFR; e.g., 
ISchiminovich et al]|2007l) and while a very weak non-linear 
trend is observed with the concentration index. The tightest cor- 
relation is with the observed NUV - r colour (Pearson correla- 
tion coefficient r = -0.85, dispersion along the y-axis cr =0.38 
dexfl while the scatter gradually increases for the stellar surface 
density (r = -0.79, cr =0.45 dex) and stellar mass (r = -0.68, 
cr =0.53 dex). These results are remarkably simil ar to the scal- 
ing relation involving the H l-to-stellar mass ratio dCortese et al.l 
1201 U ICatinella et alj|2010l; iFabello et a"l]|201 lb . suggesting that 
the dust and atomic hydrogen content of galaxies might be di- 
rectly linked. 

As discussed in the previous section, it is important to inves- 
tigate if and how the assumptions made on the dust properties of 
our sample can affect the shape and significance of the scaling re- 
lations. In Appendix A we show that only for /3 <1.5 systematic 
effects start to be significant (i.e., ~0.12 dex). However, when 
moving from y3=2 to/3 < 1 .5, the decrease in the value of the dust 
mass would be larger for massive, high stellar mass surface den- 
sity galaxies, effectively reinforcing the strong correlations be- 
tween the dust-to-stellar mass ratio and stellar mass, stellar mass 
surface density and NUV- r colours. Thus we conclude that the 
main trends shown in Fig. [3] are independent of the assumptions 
on the dust pr operties of the HRS. 

Recently, Cortese et alJ d201 li) have shown that the HI scal- 
ing relations for the whole HRS are slightly biased towards lower 
gas content, with respect to the average scaling relatio n of local 
galaxies (^g.^ lCatinella et alJl2010H Fabeno et al.l l2011l) . This is 
because, by construction, nearly half of the galaxies in the HRS 
belong to the Virgo cluster. In order to test whether this is also 
the case for the dust scaling relations, and to characterize the 
relation between dust-to-stellar mass ratio in different environ- 
ments, we have divided our sample into four different subsets: (a) 



' Concentration indexes are not corrected for inclination but this does 
not significantly affect the results presented here. 

These parameters are obtained using upper limits for the non- 
detections. 
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Fig. 3. The dust-to-stellar mass ratio as function of stellar mass (upper left), NUV-r colour (upper right), concentration index (lower 
left) and stellar mass surface density (lower right). Circles and triangles show Herschel detections and non-detections, respectively. 
The typical errors associated with our galaxies are indicated on the bottom-left corner of each panel. 



Hi-normal (i.e., Defni < 0.5, 158 galaxies) and (b) Hi-deficient 
(i.e., DefHi > 0.5, 124 galaxies) systems, (c) galaxies outside 
the Virgo cluster (144 galaxies) and (d) galaxies belonging to 
one of t he Virgo cluster clo uds (Virgo A, B, N, E and S as de- 
fined bv lGavazzi et al.ll999[ 138 galaxies). The two pairs (a)-(b) 
and (c)-(d) are mutually exclusive, while the other combinations 
are not; thus these are complementary sub-samples. The scaling 
relations for the four samples are plotted in Fig. |4] In order to 
properly quantify the difference between different environments, 
in the right most panel of Fig. |4] and in Table [T] we present the 
average trends (i.e., < log{Mdusr/Mt) >) for each subsample, 
determined by placing the non-detections to their upper-limit. 
Although for all four samples the dust-to-stellar mass ratio de- 
creases with stellar mass, colour and stellar mass surface den- 
sity, galaxies in different environments have different dust con- 
tents. In particular, for fixed stellar mass and stellar mass sur- 
face density, Virgo and Hi-deficient galaxies have, on average, 
a lower dust-to-stellar mass ratio than Hi-normal/field galaxies. 
This difference is particularly strong between Hi-normal and Hi- 
deficient galaxies (i.e., ~0.5-0.7 dex), while it is less remarkable 
(i.e., ~0.2-0.4 dex) when galaxies are separated accordingly to 
the environment they inhabit, suggesting that the atomic hydro- 



gen content is more important than the local environment in reg- 
ulating the positions of galaxies in the scaling relations. 

In order to better quantify the difference between the var- 
ious subsamples, we estimated the residuals of the Hi-deficient 
and Virgo galaxies from the mean trends observed for Hi-normal 
and galaxies outside Virgo, respectively. A simple test indi- 
cates that Hi-deficient and Hi-normal galaxies do not follow the 
same MjusilM, vs. M» and MdustlM^ vs. relations at >99.99% 
significance level. Galaxies in and outside Virgo do not follow 
the same M^„j,/M, vs. M» and MdustlM, vs. yu» relations at a 
~99.7% and ~9%% level, respectively. A similar conclusion is 
reached if, instead of performing a;^'^, we just compare the me- 
dian residuals of the two populations: while Hi-deficient and Hi- 
normal galaxies differ by ~7-8 cr, galaxies inside and outside the 
Virgo cluster show just a ~2-3 cr difference. The shift towards 
lower dust content for Hi-deficient/Virgo galaxies is also con- 
firmed by the fact that the vast majority of non-detections are 
found among gas-poor cluster galaxies. Finally, also the scat- 
ter of the Mdusi/M» vs. M, and MdusilM^ vs. ju, relations varies 
among the four sub-samples here considered: from ~0.4-0.3 dex 
to ~0.6-0.5 dex when moving from Hi-normal to Hi-deficient 
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Fig. 4. The dust-to-stellar mass ratio as function of stellar mass (top row), NUV - r colour (middle row), and stellar mass surface 
density (bottom row) for Hl-normal (first column), Hl-deficient (second column) galaxies, systems outside (third column) and 
inside (fourth column) the Virgo cluster. Circles and black triangles indicate detections and non-detections respectively. The 
average relations for Hl-normal (solid), Hl-deficient (long-dashed) galaxies, systems outside (dotted) and inside (short dashed) the 
Virgo cluster are presented in the fifth column. In the third and fourth column Hi-deficient galaxies outside the Virgo cluster and 
Hi-normal Virgo galaxies are highlighted in red and blue, respectively. 



systems and from ~0. 57-0. 44 dex to ~0.63-0.5 dex when mov- 
ing from objects outside and inside Virgo, respectively. 

It is quite easy to understand why a rough cut in environ- 
ment is less powerful in isolating dust poor systems than a cut 
in gas content. On one side, our Virgo sample still includes Hi- 
normal st ar-forming galaxies not yet aff'ected by the cluster e nvi- 
ronment (iB oselli & Gav azzi'2006VCortese & Hughes'2009). On 
the other side, we fi nd Hi-deficient galaxies also outside Virgo 
jCortese et al.ll20I lb . These two factors together reduce the dif- 
ference between the two samples. This is clearly shown in Fig.|4] 
where Hi-deficient field galaxies and Hi-normal Virgo galaxies 
are highlighted. It is also important to remember that, due to the 
lack of statistics available outside the Virgo cluster (in particu- 
lar for early-type galaxies), it is impossible to perform a more 
detailed investigation of environmental effects at the scale of 
galaxy groups and pairs. 

Contrary to the M^/„j,/M, vs. M, and Mrf„.„/M» vs. yu, rela- 
tions, the relation between M^„j,/M, vs. NUY - r does not show 



any clear variation with environment or Hi content. Although 
gas-poor Virgo systems cover a larger range of colours and are 
offset towards redder colours than Hi-normal 'field' galaxies, all 
the four sub-samples here considered apparently lie on the same 
scaling relation at a >85% significance level. This is different 
from what is observ ed in the case of the Hi scaling relations 
(ICortese etal.ll201l h. where also the Hl-fraction vs. NUY - r 
colour relation showed a variation with the environment. As we 
will see in the next sections, the independence of the M^„s,/M, 
and NIJY - r on environment is consistent with a scenario in 
which dust is removed mainly when the environment is able to 
strip the ISM directly from the optical disk, thus reducing the star 
formation and moving a galaxy along the main relation. Being 
the NUY - r colour a proxy for SSFR, our findings do not only 
confirm the recent results of da Cunha et al. ( 2010), who pointed 
out the existence of a tight relation between star formation and 
dust mass in nearby galaxies. They also suggest that this relation 
might be a very useful tool for chemical evolution models, since 
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Table 1. The average scaling relations for the four samples discussed in Sec. 3. 
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Fig. 5. The Mj„5,/M, as a function of the morphological type. 
Black filled and red empty circles are Hl-normal and HI- 
deficient galaxies, respectively. Upper limits are indicated by 
green triangles. The large empty rectangles indicate the average 
value and its error for Hl-normal (black) and Hl-deficient (red) 
galaxies in each morphological type. In case only one galaxy is 
available in a bin, only a line is shown. 

it is apparently not affected by the environment. It is important 
to remember that our colours are not corrected for internal dust 
attenuation and therefore part of the scatter at red colours could 
easily be due to the fact that we are mixing truly passive objects 
with highly obscured ones. We plan to investigate in more detail 
the relation between total SFR and dust mass in a future paper, 
once accurate dust extinction corrections are developed. 

In summary, once combined our results indicate that, at fixed 
stellar mass and stellar mass surface density, the dust content 
of nearby galaxies varies as a function of environment and Hi 
content, although the effect of the environment on the dust seems 
to be weaker than in the case of the Hi. 

Finally, since gala xy properties are also tigh tly linked to op- 
tical morphology ('e.g.. lRoberts & Havnesll994l) . it is interesting 



to investigate how the Mj„ s,/M, ratio varies as a function of mor- 
phological type. As shown in Fig.|5] the Mjusi/M, ratio rapidly 
decr eases when movin g from late- to early-type galaxies (see 
also ISkibba et al.ir201 ll) . A detailed analysis of the dust proper- 
ties of early-type galaxies in the HRS is presented by Smith et aT| 
Not surprisingly, almost all the non-detections are early- 
type galaxies. Moreover, even at fixed morphological type, we 
find that Hi-deficient objects have a significantly smaller dust 
content than Hi-normal systems, reinforcing once more one of 
the main results emerging from this section: i.e., when the Hi is 
removed from a galaxy the dust content is affected as well. This 
is even more evident when the mean M^„j,/M, ratios per bin of 
morphological type and their standard deviations (i.e., the large 
rectangles in Fig.|5]l are considered. 



4. The dust-to-HI mass ratio 

In the previous section we have shown that the M^„s,/M, ratio 
follows scaling relations very similar to the ones observed for 
the Hl-to-stellar mass ratio, suggesting that these two compo- 
nents of the ISM are tightly linked. In this section we investi- 
gate the variation of the MdustlMni ratio as a function of internal 
galaxy properties and environment following a similar approach 
as above. We focus on the separation between Hi-deficient and 
Hi-normal galaxies, since this turned out to be the ideal one to 
investigate environmental trends. 

In Fig. |6] we show the Mdusi/Mni ratio as a function of 
stellar mass, NUV - r colour and stellar mass surface density. 
Filled and empty symbols indicate Hi-normal and Hi-deficient 
galaxies, respectively. Arrows indicate both upper limits in Hi 
mass or in dust mass. Galaxies not detected at both radio and 
submm wavelengths are not shown. Contrary to what observed 
for the Mdusi/M^ ratio, the dust-to-HI ratio varies less than 2 dex. 
Almost no correlation is found with the stellar mass surface den- 
sity (r =0.1) while a weak correlation is found with stellar mass 
and NUV - r colour (r ~0.31 and dispersion along the y-axis 
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Fig. 6. The dust-to-HI mass ratio as a function of stellar mass (right panel), NUV - r colour (central panel) and stellar mass surface 
density (left panel). Filled and empty circles indicate Hl-normal and Hl-deficient galaxies, respectively. Arrows show upper limits 
in either dust or Hi mass. Typical errorbars are shown in the bottom-right corner of each panel. 



~0.37 dex). These relations become even weaker if we focus on 
Hi-normal galaxies. The MdustlMni vs. M, relation is still signif- 
icant (r ~0.28), showing that the MdustlMni ratio is an increas- 
ing function of stellar mass, a s expected from the st ellar mass vs. 
gas metallicity relation (e.g.. lTremonti et ani2004l) . Conversely, 
the correlation coefficient for the Mjust/M, vs. NUV - r relation 
drops to r ~0.07. This is mainly due to the presence of outliers 
at red colours (NUV - r >4), some of which are likely red just 
because they are highly attenuated and not because they are pas- 
sive. If we limit our analysis to galaxies with NUV - r <4 the 
correlation coefficient increases to -0.4. 

One of the most interesting results emerging from Fig. |6] 
is indeed the different behavior of Hi-deficient and Hi-normal 
galaxies. The dynamic range covered by the Mdusi/Mni ratio al- 
most doubles when Hi-deficient galaxies are included into the 
picture. At fixed stellar mass, Hi-deficient objects have system- 
atically higher MdustlMni ratios than Hi-normal galaxies. A sim- 
ilar offset is observed for the stellar mass surface density, al- 
though more pronounced at high yu», whereas in the M dust I M hi 
vs. NUV - r relation Hi-deficient galaxies appear mainly to ex- 
tend the relation already seen for Hi-normal galaxies to redder 
colours, confirming that the gas removal and the reddening of 
the colours (i.e., quenching of the star formation) are related 



qu 

dBoselli & Gavazzil2006HBosem et al.l2008HCortese & HughesI 
I2009h . This suggests that, when the gas is removed, the frac- 
tion of dust stripped from the galaxy is significantly lower than 
the Hi, thus automatically increasing the MdusilMni ratio in 
these objects. Similar results are obtained if we investigate the 
variation of the M dust I M hi ratio as a function of morphologi- 
cal type (see Fig. |7]). The Mdusi/Mni ratio is almost indepen- 
dent of morphological type, showing (in the case of Hi-normal 
galaxies) a very small decrease only towards either early-type or 
irregulars/BCDs, c onsistent with what previously observed by 
iDraine et all (12007 ) for SINGS galaxies. At fixed morphologi- 
cal type, Hi-deficient galaxies have higher M dust I M hi ratio, con- 
firming that the dust is less affected by the environment than the 
gas. 

In addition to Hi-deficient galaxies showing a higher 
M dust I M HI ratio than normal, it is interesting to note the presence 
of a few outliers having Mdust/Mni significantly lower (< lO"-') 
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Fig. 7. The M dust I M hi ratios as a function of the morphological 
type. Symbols are as in Fig.|5] 
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Fig. 8. The extended UV disk of NGC3657. SDSS-optical (left), 
GALEX-NUV (center) and //e«c/ieZ/SPlRE 250yum (right) im- 
ages are compared on the same spatial scale. 



than the average value for Hi-normal systems (~ 10 ^■'), namely: 
NGC3657, NGC4203, NGC4262, NGC4636 and UGC8045. In 
the case of UGC8045 and NGC4636 the low M dust I M hi could 
just be due to a problem with the estimate of the Hi mass. 
UGC8045 could indeed be confuse d in Hi (i.e., mor e than one 
galaxy included in the Arecibo beam. [Kent et al.l2008l) and in the 
case of NGC4636 values published in the literature vary between 
~ lO' ' Mq (Gavazzi_et_al. 2003; Kn app et al. 1978) and < 8x 10^ 
Mo dKrishna Kumar & Thonnardll983h . Conversely, we can ex- 
clude a measurement problem for NGC3657, NGC4203 and 



10 



Cortese et al.: The dust scaling relations of the Herschel Reference Survey 



NGC4262 and conclude that these systems have significantly 
higher amounts of atomic hydrogen compared to their dust con- 
tent. This is confirmed when we look in detail at the proper- 
ties of these syst ems. NGC4203 and NGC4262 have pec uliar 
Hi morphologies jKrumm et al.ll 19851: Ivan Driel et al.ll 19881) and 
several independent analyses have proposed that these galaxies 
might have recently accreted cold gas (see ICortese & HughesI 
[2009 and references therein). Although no Hi interferometric 
data are available for NGC3657, the huge ultraviolet extended 
disk revealed by GALEX images (see Fig. [8]l suggests that this 
galaxy likely contains an extended Hi reservoir. Thus, we spec- 
ulate that the unusually low M dust I M hi ratio observed in these 
systems may be due to the presence of a considerable amount 
of pristine gas not yet enriched by star formation, either because 
it has been recently accreted or because it has too low column 
density to form stars at high rate. Unfortunately, no gas metallic- 
ity estimates for the outer parts of these galaxies are available. If 
confirmed, this result would suggest the use of the MjustlMni ra- 
tio as an alternative method to identify objects that have recently 
accreted a significant amount of pristine gas. 

Finally, there is a caveat to this discussion; as discussed in 
Appendix A the relation between MdusilMni and stellar mass 
could become gradually weaker for (3 <2 (e.g., for p -1.5 
the correlation coefficient of the Mdust/Mni vs. M, relation de- 
creases to ~0.27) and it might almost disappear in the extreme 
case ofyS ~1. 



3 

T3 
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gas' \ • gas' 

Fig. 9. The dust-to-stellar mass ratio as a function of the total gas 
fractions for HRS galaxies. The solid blue line shows the predic- 
tions of a simple closed-box models (Edmunds 20 01). See the 
text for details. The five galaxies with low dust-to-Hi ratio mass 
discussed in Sec. 4 are highlighted. Symbols are as in Fig.|6l 



5. Discussion 

5.1. The origin of tine dust scaiing reiations 

We have shown that the dust-to-stellar mass ratio strongly anti- 
correlates with stellar mass, stellar mass surface density and 
colour over the whole dynamic range covered by the HRS, re- 
producing remarkably well s imilar scaling relations observed for 
the H i-to-stellar mass ratio dCatinella et al.ll20Tol : ICortese et alJ 
I2OI ih . Similarly, the dust-to-stellar mass ratio monotonically de- 
creases when moving from late- to early-type galaxies. Although 
some of these scaling relations are not new [e.g., da Cunha et al. 
(120101) has recently shown that the Mj„ „ /M, ratio strongly anti- 
correlates with SSFR], this is the first time that we can accurately 
quantify them for a volume-limited sample of galaxies spanning 
all morphologies and for which Hi information is available. 

Once the prediction of dust formation an d evolution mod- 
els are taken into account (e.g., iDwekl FT998t lEdmun dsl l200H 
ICalura et al.ll2008l) . it is quite easy to understand where these 
scaUng relations are coming from. As nicely described by 
Ida Cunha et alj (1201 Ol) . the trends shown in Fig. [3] can be seen 
as the result of th e variation of SSFR with stellar mass (e.g., 
iBoselli et aDl200ll: ISchiminovich et al.ll2007l) . Low mass, late- 
type systems are characterized by high SSFR and gas fraction. 
The gas can sustain the star formation activity and, consequently, 
a large fraction of dust is formed, likely exceeding the amount 
of dust grains destroyed in the ISM. When we move to higher 
stellar mass and earlier morphological types, the SSFR and gas 
fraction start decreasing and, as a result, the amount of dust pro- 
duced will no longer be able to overcome those destroyed in the 
meantime. This would naturally explain why the dust-to-stellar 
mass ratio decreases with stellar mass, specific star formation 
rate, stellar mass surface density and when moving from irreg- 
ulars to elliptical galaxies. In other words, what we are looking 
at is the change of the amount of baryonic mass stored in the 
ISM as a function of internal galaxy properties, explaining why 
Hi and dust follow the same scaling relations. 



To test this scenario qualitatively, w e use the simple analyt- 
ical approach developed by .Edmundsl (1200 lb in the case of a 
closed-box. We use his Eq. 5 and 1 1 to determine how the dust- 
to-stellar mass ratio evolves as a function of the gas fraction. The 
model has six free parameters: namely, the efficiency of dust 
condensation from heavy elements made in stellar winds from 
massive stars and supernovae (xi) and from Asymptotic Giant 
Branch stars {X2), the effective yield (p), the fraction of mass 
of the ISM locked into stars (a), the fraction of the ISM where 
mantles can grow (e) and the fraction of dust d estroyed by star 
formation (6). Here we assume vi - 0.2 ( Mor gan & Edmundi 
2003i) . Y-, =0.5 dZhukoyska et al.[l 2008). p = 0.012 an d a =0.80 
(i.e., consistent with a lChabried l2003 IMF), 6 =0.3 (lEdmundsl 
2001') and e =0.5. In order to estimate the total amount of gas 
of our galaxies from the Hi masses, we assume a molecula r-to- 
atomic hydrogen gas ratio of 0.38 dSaintonge et al.ll201 ll) and 
include the contribution of helium and heavy elements (i.e., a 
factor 1.36). The overall trend is not significantly affected if we 
account only for the contribution from atomic gas in the estimate 
of the gas fraction. We note that the goal of this exercise is just to 
establish whether the predictions of a simple closed-box model 
are roughly consistent with the shape of the trends we find. There 
are many free parameters and assumptions in the modeling, as 
well as a large uncertainty in the absolute value of our dust 
masses, thus looking for an exact fit does not seem meaningful 
at this stage. In addition, this is just a simple close-box model, 
whereas the evolutionary history of some galaxies in our sam- 
ple (in particular early-type systems) might have been signifi- 
cantly different. The prediction of this simple model are shown 
in Fig.|9] As expected, the Mdust/M, ratio decreases with the gas 
fraction confirming that the scaling relations we found are just a 
consequence of the tight link between dust and the cold gas in 
the ISM. Interestingly, the five galaxies with unusually low dust- 
to-Hi ratio, discussed in § 4, are also outliers in Fig.|9] Thus, for 
these objects, a closed-box approximation is likely inadequate. 
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Our model is also unable to reproduce the very dust-rich sys- 
tems in our sample. A similar tension between observations and 
models was recently found by Dunne et al. (2011) studying the 
properties of galaxies in the H-ATLAS survey. By comparing 
more detailed models predictions with observations they show 
how theoretical models are not able to reproduce galaxies with 
\og{MdustlMt) > -2.5. The origin of this disagreement is still 
unclear, but we note that it may partially be solved if /3 is lower 
than 2. We plan to investigate further this issue in the future when 
more accurate dust estimates will become possible. 

5.2. Dust stripping 

The other important result emerging from our analysis is that at 
fixed stellar mass, stellar mass surface density and morphology, 
Hi-deficient/Virgo cluster galaxies have systematically lower 
Mdusi/M» and higher M,;„5,/M/// ratios than Hi-normal/field sys- 
tems. Once combined, these two independent results provide a 
consistent picture of the effects of the environment on the dust 
content of cluster galaxies. On one side, the reduction in the 
Mch(st/Mt ratio is a clear evidence for dust removal in cluster 
galaxies. This is reinforced by Fig. [10] (left), where we show 
that the Mj„5,/M» r atio strongly decrea ses as a function of Hi 
deficiency (see also ICorbelli et al.ll201 ih . On the other side, the 
higher M dust I M hi ratio indicates that the effect of the environ- 
ment on the dust content is significantly less dramatic compared 
to the amount of Hi stripped from the disk and dispersed into 
the i ntra-cluster medium. As akeady suggested bv lCortese et al.l 
(l201 0b). this could be easily explained if the dust disk is sig- 
nificantly less extended than the gas disk (iThomas et al. 2004*) 
or, in other words, if the dust-to-gas ratio decreases mono- 
toni cally with increasing distance from the cen t er of galax- 
ies ('Bendo et al."2010'; 'Mun oz-Mateos et al.ll2009t iPohlen et all 
12010; Magrini et al. 2011). 

We test this interpretation in the right panel of Fig. (TO] 
where we compare the strong correlation between MdusilMt and 
Mni/Mt found in our data, with the predictions of a very simple 
toy model. We assume that the stripping of the ISM is outside- 
in, with the outer parts being affected before the galaxy center, as 
expected in a ram-pressure str i pping scenari o (Koocmann et al. 
2006HBoselli & Gavazzil2006tBoselli et al.ll2006: .Cortese et al.i 



201 In . Once the stripping starts, all the components of the ISM 
(i.e., both gas and dust) are completely removed from the disk 
up to the truncation radius. The dust mass surface density profile 
is exponential while we consider two different surface density 
distributions for the Hi. In the first one (solid line in Fig. [TOt . 
the Hi profile is exponential and the sca le length of the Hi dis k is 
just 1 .5 times the one of the dust disk dThomas et al. I I2OOI . In 
the second on e (dashed line in Fig.fTOll. the Hi profile is constant 
(iBidel etal.l [2010) and extends up to 8 times the scalelength of 
the dust disk. This corresponds to a Hi radius 1 .5 times the op- 
tical/dust radius assuming that the dust follows the stellar distri- 
bution and that the ratio betwe en optical scale-length and optical 
radius is -5.3 (iGavazzi et al. I I2OOQ) . Both models are normal- 
ized so that the unperturbed galaxy has the typical log(Mrf„5, /M.) 

of Hi-normal galaxies ( 2.46 dex) and Mdusi/Mni =100. We 

note that these two cases represent two extremes in the differ- 
ence between the dust and Hi distribution in galaxies. In Fig.fTOl 
(right) we show the evolution of the Mjusi/Mt and Mhi/M, pre- 
dicted by the two models for different values of stripping radius. 
As expected, the range covered by these models reproduces the 
difference in Mjusi/Mt and MhiIM^ between Hi-normal and Hi- 
deficient galaxies. This supports our simple interpretation that 
the dust is affected by the environment in a less dramatic way 



than the Hi just because, like the molecular gas, it is more cen- 
trally concentrated than the atomic gas component. These re- 
sults confirm once more dust stripping by environmental effects 
as an important mechanism for injecting dust grains into the 
intra-cluster medium , thus contributing to its metal enrichment 
(iPopescu et al.ll2000l) . 

The clear evidence of dust stripping emerging from our anal- 
ysis leaves open the possibility that also the molecular hydro- 
gen component in galaxies could be significantly perturbed in 
the center of cl usters of galaxies. However, this is sue is still 
highly debated (iBoselli et al.ll2002l: iFumagalli et"al] 2009) and, 
as clearly shown by our analysis, only a detailed comparison of 
molecular hydrogen properties for galaxies in different environ- 
ment but with the same morphological type and stellar mass, will 
eventually allow us to quantify if and how molecules are directly 
removed from the star forming disks of infalling spirals. 

6. Summary 

In this paper, we have combined i/ersc/ie //SPIRE observations 
with optical. Hi and UV data to investigate the relation between 
dust content, galaxy properties and environment for the Herschel 
Reference Survey, a volume- magnitude- limited sample of -300 
galaxies. Our main results are as follows: 

- We find that the dust-to-stellar mass ratio strongly anti- 
correlates with stellar mass, stellar mass surface density 
and NUV - r color across the whole range of parameters 
covered by our sample. These relations are followed by 
all galaxies, regardless of environment or gas content, but 
Virgo/Hi-deficient galaxies show a systematically lower 
dust-to-stellar mass at fixed stellar mass and morphological 
type. 

- Significantly weaker correlations are found between the 
dust-to-Hi mass ratio and internal galaxy properties, with 
the Mdusr/Mffj ratio mildly increasing with stellar mass, 
as expected from the mass-metallicity relation. Gas-poor 
cluster galaxies have systematically higher dust content per 
unit of Hi mass than Hi-normal systems. 

- We show that the differences observed between Hi-deficient 
and Hi-normal galaxies provide strong evidence for dust 
stripping in cluster galaxies. However, we demonstrate that 
the fraction of dust removed from the disk is significantly 
lower than the Hi, and that this is likely just a consequence 
of the fact that the Hi disk is much more extended than 
the dust disk or, in other words, that the dust-to-atomic 
hydrogen ratio quickly declines in the outer regions. 

- Finally, we compare our results with the prediction of sim- 
ple models of dust formation and evolution showing that the 
trends here presented are consistent with a simple picture in 
which the amount of dust in galaxies is regulated by the star 
formation activity and cold gas content. 

Although this work represents just the first step in the under- 
standing of dust properties in galaxies, our investigation high- 
lights the power of the HRS as an ideal local sample for galaxy 
evolution studies. The natural extension of the present analysis 
will be the study of the relation between dust properties, molec- 
ular hydrogen content and gas metallicity. Only after molecules 
and heavy elements in the ISM are included into the picture and 
more accurate dust mass estimates are available, will it be re- 
ally possible to unveil the role of dust on galaxy evolution not 
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Fig. 10. Left: The relation between dust-to-stellar mass and HI deficiency. Right: The relation between dust-to-stellar mass and HI- 
to-stellar mass ratio. The dotted lines indicate constant Hl-to-dust ratios of (from left to right): 30-100-300. The green lines indicate 
two different models of dust stripping assuming a constant (dashed) or exponential (solid) Hi surface density profile. See text for 
details. Symbols are as in Fig.|6] 



only providing strong constraints for theoretical models but also 
a reference for high-redshift studies, i.e., the primary goal of a 
survey such as the HRS. 
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Appendix A: Testing the dust mass estimate based 
on SPIRE flux densities only. 

In this Appendix we investigate the reHabihty of the method de- 
scribed in § l2.2l to estimate dust masses. We proceed as follows. 
We compare our estimates with different methods all based on 
roughly the same assumptions on the properties of the underly- 
ing dust population (i.e., the value of fi and its variation with fre- 
quency). This is aimed at validating the methodology presented 
here. Then, we use our method for different values of fi to study 
how the assumption fi-2 made in this paper could affect our re- 
sults. In both cases we are interested in the average difference 
and standard deviation between the various methods and in the 
presence of any possible systematic trend with stellar mass, stel- 
lar mass surface density and NUV - r colour 

The results of the first step are shown in Fig. I A. II In each 
row we show the difference between the logarithm of the dust 
mass obtained from two different methods as a function of stellar 
mass, NUV - r colour and stellar mass surface density. In the 
right most panel we also present the histogram of the distribution 
of the difference, its average value and standard deviation. 

First of all, we com pare our results with the ones obtained 
bv Davie s et alj (1201 2h for the 55 galaxies in common with our 
work (top row in Fig. lA.lb . These are obtained using a single 
modified black-body SED fitting of the five PACS-nSPIRE flux 
densities at 100, 160, 250, 350 and 500 /im and assuming the 
same dust emissivity as in this paper. In our analysis we only 
considered the points for which the SED fitting routine gives 
a. < 1, but in Fig. lA.ll we also show as empty circles the 
galaxies with > ^- We do not find any systematic difference 
in the average dust mass estimate and the scatter between the 
two methods is of the order of ~0. 1 dex, not significantly larger 
than the uncertainty of ~0.08 dex in the dus t mass estimate from 
the SED fitting technique ( Davies et al.ll201 2). Small systematic 
trends in the dust mass estimate are found as a function of stellar 
mass and stellar mass surface density. However, these are just 
marginal and within the error assumed throughout the paper 

Secondly, we consider the case of SED fitting of SPIRE- 
only points, following the same technique described above. 
Interestingly, we find a perfect agreement between the two meth- 
ods and we could in theory just use the SED fitting of SPIRE 
fluxes only. The main reason why we decided not to do so 
is that the best fi tting SED is signifi cantly different from the 
one obtained by iDavies et al.l (1201 2h when combining PACS 



14 



Cortese et al.: The dust scaling relations of the Herschel Reference Survey 



Q 



Da 



Q 



T3 oa 



O _ 
< 3 

Q 

I 



pa 



J' 



Da 
Da 



^ -1 



I ' ' ' ' I ' ' ' ' I 



o 



I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I 



_|_L_1_1_1_I_1_1_1_1_L 



I ' ' ' I ' ' ' I ' +' ' I ' ' ' ' I ' ' ' ' I 



off. 



I I I I I I I I I 



-I— H- 




_l_i_i_i_L 



J I 1 1 1 1 1 1 1 1 1 1 1 1 



I I I I I I I I 1 .1 1 1 1 1 1 1 1 1 1 1 1 1 1 . 



H T I 1 1 1 1 1 1 1 1 1 1 1 r 







8 9 10 11 
log(M.) [MJ 



2 4 6 
NUV-r [mag] 



7 8 9 IC 
log(M.) [M^ kpc-2] 







Fig.A.l. The logarithmic difference between the dust mass estimates presented in this work and those obtained from modified 
black-body SED fitting of PACS +SPIRE (top row) and SPIRE-only (second row) data, using the SED libraiy of Draine & Li (2007, 
third row) and using just the 350//m luminosity without any colour information (bottom row). For each row, the right panel shows 
the histogram of the difference, its average value and standard deviation of the population. The empty circles in the top two rows 
indicate points for which the SED fitting routine provided ax^ > 1 • 



and SPIRE data since no constraint is provided below 250//m. 
Therefore these are likely not a good representation of the far- 
infrared/submm emission for our galaxies and we preferred not 
to use them. Moreover, it is sometimes difficult to trust the un- 
certainty determined from the fitting technique since only based 
on the three SPIRE data. 

Thirdly, we compare the estimate obtained from our method 
assuming a modified black-body SED and the one obtained from 
the same method but assuming the Draine & Li (2007) SED 
(DL07) as described in § 12.21 As shown in the third row of 
Fig. lA.ll we find again a very good agreement between the two 
methods. This is due to the fact that, in th e wavelength range 
covered by SPIRE, the SED predicted by the Draine & Li (2007) 
models is very similar t o a single ternperatu re modified black- 
body with /3=2 (see also lMagrini et al.ll201 ih . 

Finally, we consider the very simple case where only a 
submm flux density is available and no colour can be used (bot- 
tom row in Fig. lA.ll ). In details, we just assumed the same 250- 
to-500 yum ratio for all our galaxies (i.e., the mean value for the 



detections in our sample /og(/25o//50o) =0.75). This is done to 
investigate the role played by the colour term in our method and 
to warn how this very rough technique could affect scaling re- 
lations studies, in particular at high-redshift. As expected, this 
is by far the worse method to estimate dust masses showing a 
scatter of ~0.2 dex. More importantly, we find systematic trends 
with stellar mass and stellar mass surface densities showing that, 
when compared with our method, this technique underestimates 
dust masses for low-mass galaxies and overestimate those of 
massive systems. The reason for these systematic trends is due 
to the fact that the 250-to-500 jum flux density ratio strongly cor- 
relates with stellar mass and s tellar mass surface density (see 
Fig.lA:2land lBoselli et al.l2012h . Thus not accounting for the dif- 
ferent colours in our galaxies introduces systematic effects that 
can be seen in Fig. lA.3l where we compare some of the scaling 
relations discussed in this paper (orange dashed regions) with 
those obtained using this method. While, as expected, the rela- 
tion with NUV - r colour is not affected, the trends between 
dust-to-stellar mass ratio, stellar mass and stellar mass surface 
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Fig. A.2. The 250-to-500 //m flux density ratio as a function of 
stellar mass (left), NUV - r colour (center) and stellar mass sur- 
face density (right). 



density change. We thus suggest caution in using a technique 
based on only a single luminosity for a detailed quantification of 
dust properties in galaxies. 

In summary, these first tests confirm that our technique does 
not introduce any significant systematic effect, once the proper- 
ties of the underlying dust component are fixed. 

The next step is to see what happens if we change the 
value of /3, keeping fixed the value of dust opacity at 350yL(m. 
In the following, we consider modified black-bodies with 1< 
/? <2.5 which appear the most favoured range of values ac- 
cording to theoretical models and observatio ns ( Finkbeiner et al] 
1999'; Duoac et al. 2003: Gordon et al.l2010 ':?aradi s et al.l201oF 
Planck Collab oration et al. 2()ljj). It is important to note that, 
given the strong degeneracy between /? and dust temperature 
we had to modify the range of dust temperatures investigated 
in order to recover the range of submm colours observed in the 
HRS. In detail, while fory6=2 and 2.5 we varied the tempera- 
ture between 5 and 55 K, for /5-1.5 and 1 we had to vary the 
temperature between 15 and 230 K. Moreover, the case y3=l is 
not able to reproduce the bluest colours observed in our sample 
(i.e., /250//500 > 7.5, 37 galaxies in total) suggesting that this 
case is not a good representation of the dust SED for at least a 
fraction of our sample (see Fig. IA.4b . The results of our inves- 
tigations are summarized in Fig. lA.5l As expected, the absolute 
value of the dust mass varies with /3 by almost a factor 10. Values 
of y6 lower than 2 give dust masses smaller than the one used in 
this work, and viceversa. However, once this offset is taken into 
account, the intrinsic scatter in the dust mass estimate stays be- 
low ~0.1 dex for 1.5< /3 <2.5 and goes up to just 0.17 dex for 
/?=1 (for this case only those galaxies with /250//500 < 7.5 are 
included), although with a significant tail of outliers. Moreover, 
very mild systematic trends are observed as a function of stel- 
lar mass, colour and stellar mass surface density but these could 
become important only for /3 ~ 1 . As discussed in the text, this 
would reinforce the trend observed between the dust-to-stellar 
mass ratio but weaken the one between dust-to-Hi mass ratio 
and stellar mass. 

In conclusion, we can confidently state that the method pre- 
sented here to estimate dust masses is a reliable solution when 
accurate SED fitting is not possible. The typical scatter with re- 
spect to our methods is of the order of ~0.2 dex. Although dif- 
ferent values of /3 would imply different values of dust masses, 
the scaling relations discussed in this paper should remain valid 
atleastforl.5<yS<2.5. 
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Fig.A.3. The dust-to-stellar mass ratio as a function of stellar 
mass (left), NUV - r colour (center) and stellar mass surface 
density (right), for dust masses estimated from the 350/vm lu- 
minosities without taking into account any colour/temperature 
variation in the sample. Symbols are as in Fig. 3. The orange 
lines and dashed regions show the average scaling relations and 
+ lcr area presented in Fig. 3. 
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Fig. A.4. The relation between the ratio MdustlifaoD^) and the 
250//m-to-500//m flux density ratio for a modified black-body 
SED withy8=2.5 (blue), 2 (black), 1.5 (green) and 1 (magenta). 

Appendix B: Recipes to estimate dust masses with 
SPIRE data only 

In this Appendix we provide the polynomial fit to the relations 
between the ratio M^uisi/ihsoD^) in units of kg W"' Hz and 
the 250jum-to-500yL(m flux density ratio which we used to esti- 
mate dust masses. We provide two different sets of coefficients 
obtained by convolving modified black-bodies characterized by 
different values of /3 (i.e., 1, 1.5, 2 and 2.5) with the SPIRE 
RSRF for extended and point-like sources (Table IBTI ). Once the 
value MdiistlihsoD^) has been computed, the dust mass can be 
obtained from the following equation 

We note that the difference degrees used to fit fi >2 and /3 <2 
are due to the different range of temperatures adopted for each 
case. For /3 <2 we had to use larger temperature range (15 to 
230 K instead of 5 to 55 K) and, as a consequence, an higher or- 
der polynomial to fit the relation between the 250//m-to-500jum 
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Fig.A.5. The logarithmic difference between the dust mass estimates presented in this work and the one obtained assuming a 
modified black-body with 13-2.5 (top row), 1 .5 (middle row) and 1 (bottom row). For each row, the right panel shows the histogram 
of the difference, its average value and standard deviation of the population. 



and Mjust/ifssoD^) ratios. The fits are optimized for the range of 
colours covered by our sample, i.e., 3< /250//500 <11- 
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Table B.l. Relations to determine the MdustlihioD^) from the 250/zm-to-500//m flux density ratio for extended sources and point 
sources. 
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Point Sources 
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log(/250//500) 
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-284.971 
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lOg(/250//500) 


16.171 
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-110.479 
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log(/250//500) 


16.829 
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-0.116 
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log(/250//500) 


17.053 
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Extended Sources 
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